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' In this paper we consider a system consisting of a two-level atom in an excited state inter- 

bJ)' acting with two modes of a radiation field prepared initially in ^-photon coherent states. This 

system is described by two-mode multiphoton (, i.e., ki, k^) Jaynes-Cummings model (JCM). 
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For this system we investigate the occurrence of the revival-collapse phenomenon (RCP) in 
the evolution of the single-mode, two-mode, sum and difference quadrature squeezing. We 
show that there is a class of states for which all these types of squeezing exhibit RCP similar 
to that involved in the corresponding atomic inversion. Also we show numerically that the 
single-mode squeezing of the first mode for (k\,k2) = (3,1) provides RCP similar to that 
of the atomic inversion of the case (kx,k2) = (1, 1), however, sum and difference squeezing 
give partial information on that case. Moreover, we show that single-mode, two-mode and 
sum squeezing for the case (fci,^) = (2,2) provide information on the atomic inversion of 
the single-mode two-photon JCM. We derive the rescaled squeezing factors giving accurate 
information on the atomic inversion for all cases. The consequences of these results are that 
the homodyne and heterodyne detectors can be used to detect the RCP for the two-mode 
JCM. 



> 

On 
00 



oo 
o 

O ■ PACS numbers: 42.50.Dv, 32.80.-t, 42.50.-p 



X 



I. INTRODUCTION 

The simplest model in quantum optics is the Jaynes-Cummings model (JCM), in which a 
radiation field interacts with a single two-level atom [l[. This system has become experimentally 
realizable with the Rydberg atoms in high-Q microwave cavities (, e.g., see In the framework of 
the rotating wave approximation many of interesting effects have been reported for JCM. The most 
important phenomenon is the revival-collapse phenomenon (RCP), which occurs in the evolution 
of the atomic inversion (a z (T)), i.e. instead of displaying steady Rabi oscillations in the case of a 
classical field coupled to the atom [3j, there is an initial collapse of these oscillations followed by 



regular revivals that slowly become broader and eventually overlap 



a. 



This indicates that RCP 



is a pure quantum mechanical effect [5j. For more details about the RCP of the JCM the reader 



can consult, e.g. [4]. 

The JCM has been extended to include multimode fields, e.g. [f], B, multilevel atoms 
and multiatom interactions [ldl. Two-mode JCM (TJCM) has taken a considerable interest and 

n □ nu 

studied from different points of view, e.g. [6|, |8|, llll, 1121]. RCP of the TJCM is rather complicated 
compared to that of the single-mode JCM in the sense that the revival series is compact and each 
revival is followed by secondary revival. Furthermore, for strong intensities the locations of the 
revival patterns in the time interaction "domain" are independent of the intensities (see Fig. 1). 

Ill ], however, an investigation for the occurrence 



Such behaviour has been partially explained in 
of the secondary revivals is given in Q]- 

Quite recently a new technique is developed for discussing how within the single-mode multipho- 
ton JCM the RCP of the atomic inversion of the standard (, i.e. single-photon) JCM is manifested 
in the evolution of the quadrature squeezing of the field [13[] . Two approaches have been adopted 
for such analysis, namely, natural and numerical-simulation approaches. For natural approach it 
has been shown that there is a class of states whose squeezing factors can directly include informa- 
tion on the corresponding atomic inversion. However, the numerical-simulation approach has been 
given to show that the evolution of the quadrature squeezing of the three-photon JCM reflects the 
RCP involved in the (a z (T)) of the single-photon JCM for the same initial field state. Moreover, 
we have deduced a general form for the higher-oder squeezing factor, which can give information 
on the atomic inversion of the single-photon JC M [l^ and two-photon JCM [lj|. 

In this paper we apply the technique given in [13|] to TJCM for investigating the occurrence of 
the RCP in the quadrature squeezing and how can be connected with the atomic inversion. For 
convenience we assume that the radiation fields are initially prepared in /-photon coherent states 
[if] , [it], [rJ and the atom is in an excited atomic state. Needless to say that the situation for the 
TJCM is more complicated than that of the single-mode JCM. For instance, there are different types 
of quadrature squeezing such as single-mode, two-mode, sum and difference squeezing. Moreover, 
the strong entanglement between the two bosonic systems over the atomic system making the 
investigation is rather complicated. In spite of these difficulties we have obtained many interesting 
results, e.g. for all types of squeezing there is a class of states for which squeezing factors can 
directly give the corresponding atomic inversion. Additionally, using numerical technique we have 
shown that when k\ + k% = 4 (cf. (QJJ) the F-quadrature squeezing factor of the particular types 
can provide RCP similar to that exhibited in the evolution of the atomic inversion of the standard 
(, i.e., k\ = &2 = 1) TJCM or single-mode two-photon JCM based on the values of kj. We 
have to stress that the nonclassical squeezing for TJCM has been studied by several authors, 



e.g., see [191, and it will not be considered in the present paper. Finally, the results given here 

rfnn 

and in [13j, LL4|, ll5f] show that the RCP occurred in (a z (T)} can be detected using techniques 
similar to those used for quadrature squeezing, e.g. homodyne detector j^J], nonlinear homodyne 

22l | . It is worth mentioning that in cavity QED, 



detector [21[ and multiport homodyne detector 



the homodyne detector technique has been applied to the single Rydberg atom and one-photon 
field for studying the field phase evolution of the regular JCM {23]. Quite recently similar s etup 



is given for induced measurement and quantum computation with atoms in optical cavities 

n 



2J. 



Moreover, the progress in both of the trapped ions [25] and micromaser [26] is promising to produce 
the phenomena discussed in the paper. We conclude this part by drawing the attention to that 
the first experimentally observed squeezed states are of the two- mode type j^ . 

The paper is prepared in the following order. In section 2 we give the basic relations and 
equations including the model and the definition of squeezing. In sections 3-5 we investigate 
single-mode squeezing, two-mode squeezing and sum-difference squeezing, respectively. In section 
6 the main conclusions are summarized. 

II. BASIC RELATIONS AND EQUATIONS 

In this section we give the basic relations and equations, which will be used throughout the 
paper. Precisely, we write down the Hamiltonian of the system under consideration, its wave 
function and the definition of quadrature squeezing. 

n 

The Hamiltonian controlling the TJCM in the rotating wave approximation is [121 ]: 

/./ i ^ A + A 1 A / « k-i ~ki * "ffci~ iko \ / \ 

— = LOia[ai + W2« 2 a 2 + 2 LOa<Tz + 2 + + 1 2 <T ~'' ' ' 

where a± and a z are the Pauli spin operators; Uj, (j = 1, 2) and uj a are the frequencies of the cavity 
modes aj and the atomic frequency, respectively; g is the atom-field coupling constant and kj is 
the transition parameter of the jth mode. The derivation of the Hamiltonian (pTJ) from the first 
principle is given in [281 ] . 



We restrict the investigation to the exact resonance case k±ui + k%ux = Wa- For evaluating the 
dynamical state of (PQ) we define two operators F\ and F2 as 

Fi = u)ia\ai + cj2«2"2 + ^w a <r z , F 2 = g(a kl a2 2 a + + a\ kl a^ 2 <r_ ) . (2) 

It is easy to prove that F\ and F2 are constants of motion. This fact leads to that the evolution of 
the mean-photon number of the modes and the atomic inversion of the system include information 
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on each other. In the interaction picture the unitary evolution operator of the Hamiltonian ([I]) 
takes the form 

Ui(T,0) = exp(-ifF 2 ) 



cos(TD) - i^^F 2 



(3) 

-2- 

where 

T = gt, D 2 = a^a k Ja\ kl a\ k2 d + a„ + a\ kl a\ k2 a k ^a k ?d„a+. (4) 

For the sake of generalization we consider the jth mode is initially prepared in the /-photon 
16. iirl. \\\ having the form 



coherent states 



00 1 a n 

1^(0)) = E^IW. ^' ) = exp(--|a J | 2 )-^, (5) 



v'n! 

ri=0 



where lj are parameters their values will be specified in the text. Also throughout the paper 
we consider a,- are real. States (EI) can be obtained from Ith harmonic generation using Brandt- 

n 

Greenberg operators [18]. We proceed by considering that the atom is initially in the excited state 
|+). Therefore, the total initial state of the system is 

l*(0)> = |^i(0)>(g)|^(0)>(g)|+>. (6) 



From ([3]) and ([6]) the dynamical state vector of the system can be evaluated as 

oo 

|*(T)) = C ntm [cos(TA njm )\ + J 1 nJ 2 rn)-ism(TA ntm )\-J 1 n + ki,l2'rn + k2)}, (7) 

n,m=0 

where |— ) denotes ground atomic state, C n ^ m = Cn^Cm^ and 



_ / (hn + fci)!(/ 2 m + k 2 )l . . 

n ' m "V (hn)Khm)\ ■ W 

The atomic inversion associated with © is 

oo 

(a z (T)) = J2clm^ S (2TA n , m ). (9) 

n,m 
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FIG. 1: The atomic inversion (& Z (T)) against the scaled time T when the optical cavity modes are initially 
prepared in the coherent states for (k\, fe, ot\, 012) = (1, 1, 5, 5). 

To investigate the evolution of quadrature squeezing we evaluate the general form for the dif- 
ferent moments of the a] and aj for ((TJ) when 1% = I2 = 1 as 



1 / / 1 / / OO 

(«l Sl ( r )ai 2 ( r )«2 S3 ( T ) 6 2 4 ( r )) = E C f n+ a / 1 ,«H-4 C 'n+4,m+«i 

n,m=0 



cos(TA n+s / jm+s / ) cos(TA, 



n+s 9 ,m+s. 



V(»+4)!(n+4)!(m+s^)!(m+4)! 



(10) 



J-oWTA WWTA ^ x /(n+fc 1 +s' 1 )!(n+fc 1 +4)!(m+fc 2 +4)!(m+A ;2 +4)! 

"I" Sm^l ii. n+s / ijm+ ^ j SlIHI A n+S ^ m+S ^ (n+fei)!(m+fe 2 )! 



where Spj = 1,2,3,4 are positive integers. Also we define two quadratures X and Y, which 
denote the real (electric) and imaginary (magnetic) parts of the radiation field. Assuming that 
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these quadratures satisfy the following commutation rule: 

M = y, (11) 

where d may be c-number or operator. The uncertainty relation related to the commutation rule 
CEH) is 

((AA) 2 )((AY) 2 >>MJ! (12) 

where ((AX) 2 ) = (X 2 ) — (X) 2 and similar form can be given for ((Ay) 2 ). The system is said to 
be squeezed in the A-quadrature if 

S(T) = 2((AX(T)) 2 )-±\(d)\<0. (13) 

The equality sign in (|13p holds for minimum-uncertainty states. Similar definition can be given for 
the Y-quadrature (defining a Q-factor). As we mentioned in the Introduction we study the evolu- 
tion of four types of quadrature squeezing: single-mode, two-mode, sum and difference squeezing. 
The object of such study is to follow the possible occurrence of the RCP in the evolution of the 
squeezing factors and the conditions required for such occurrence. Also we try to find which type 
of squeezing factors can fit well information on the evolution of the atomic inversion. These issues 
will be discussed in the following sections. 

III. SINGLE-MODE SQUEEZING 

In this section we study the occurrence of the RCP in the evolution of the single-mode squeezing 
factors for TJCM. The single-mode squeezing factors Sj(t) and Qj(t) for the jth-mode when the 
quadratures X and Y are defined in the standard form, can be expressed as 

S^T) = (a](T)aj(T)) +Re(d 2 (T)) - 2(Re(a,(T))) 

(14) 

Q 3 (T) = (aj^a^T)) -Re(a 2 (T)) -2(lm(a J (T))) 2 . 
In the following parts we investigate the natural and numerical approaches in a greater details. 

A. Natural approach 



In this part and throughout the paper the natural approach is given for the standard TJCM. 
This approach is based on the fact that F\ is a constant of motion and hence the quantities 
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FIG. 2: The single-mode squeezing factors against the scaled time T when the modes are initially prepared 
in the coherent light with (ai,a2) = (5,5) and for (kijkz) = (3,1) (a), (1,3) (b) and (2,2) (c). 

(a z (T)), (a\(T)ai(T)) and (o2(T)ct2(T)) can carry information on each other. Thus the squeezing 
factors (j!4j) can give information on (a z (T)) when 



(a i (T))=0, (a J 2 (T)) = 



(15) 



simultaneously. This situation can be established when the jth mode is initially prepared in three- 



photon states 2(]], four-photon states 30(] and so on. Also for particular values of the parameter 
I the /-coherent state © can fulfill conditions (|15p . For such type of initial states one can easily 
prove that 

(a z (T)) = 2<ot(0)of(0)) + 1 - 2Sj(T). (16) 
Expression (|16|) shows that the atomic inversion can be readout from the quadrature squeezing. 



B. Numerical simulation 



In this part and throughout the paper numerical-simulation approach is applied to TJCM when 
k\ + ^2 > 2 and the modes are initially prepared in coherent light, i.e. l\ = I2 = 1 in Now 
the object here is to discuss the possibility of obtaining RCP in the evolution of the single-mode 
squeezing factors similar to that of the atomic inversion of the standard TJCM, i.e. (<5" 2 (T))fc 1= fc 2= i. 
The procedures related to this technique are given in [7| and we briefly explain them for the first 
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mode. From ()14[) RCP may occur in Sj(T) (or Qj(T)) only when the evolution of the Re(a :) '(T)} (or 
Im(dj(T)}) are close to zero (, i.e. steady state) since these quantities are squared. Thus when the 
probability amplitudes C n ^ m are real Im(aj(T)) = and consequently the RCP can likely occur in 
the evolution of Qj(T). Additionally, when k\ + &2 > 2, (aj(T)aj(T)) exhibits chaotic behaviour. 
From numerical data for this case we can consider {a]{T)aj{T)) ~ (at(O)aj(O)). This means that 
the occurrence of the RCP (if it is so) in Qj(T) is related to the quantity Re(a 2 (T)). Consequently, 
we compare the form of Re(a|(T)) with that of (a z (T))) cl= i i . 2= i. Now we give a closer look at 
(a 2 (T)), which from (jlOp has the form: 



(aj(T)) = aj Z P(n)P(m) cos(TA n+2im ) cos(TA n , m ) 

n,m=0 



(17) 



( n+fcl+ l)( w +fc 2 +2) a ^frp A ,„ ^ s i n (TA nr 



(n+l)(n 



±ff^sin(TA 



n+2,m ) 



where P(n) 



(1)\2 



is the photon- number distribution for the coherent light. In the strong- 



intensity regime, i.e. fij = (dj(O)aj(O)) = |ay| 2 >> 1, and finite values of the transition parameters 
kj we can apply the harmonic approximation technique j^, 11]. This technique is based on the 
fact that the photon-number distribution of the coherent light is Poissonian with a sharp peak at 
n = n and hence the terms which contribute effectively to the summation in (|17|) are those for 
which n ~ n. As a result of this fact the square root in the second line of (|17p tends to unity and 
(fT?]) reduces to 



(a 2 (T))~ni J2 P(n)P(m)cos[T(A n+2 , m -A n)m )}. (18) 

n,m=0 

Regardless of the prefactor fi\ in (fT8"|). the comparison between ([9]) of the case (kj,lj) = (1,1) 
and (|18p leads to that both expressions exhibit quite similar dynamical behaviour only when the 
arguments of cosines in the two expressions are comparable. Therefore, we seek the proportionality 
factor fix, say, which can be evaluated from the following expression 

An+2,m A n m . > 

III = ; (19) 

2^(n + l)(m + l) 
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FIG. 3: The modified squeezing factor for the cases (fci, ki) — (3, 1) (a) and (fci, k-z) = (2, 2) (b) when the 
modes are initially prepared in the coherent light with ot\ = a.2 — 5. 



Expression (|19|) can be re-expressed as 



Ml 



(m+fc 2 )K ra +fcl 



kiV- f 

!) ! \2j(n+2)(n 



(2n+3)fci+fcf 



(n+l)!(m+l)! \ 2 ^/(n+2)(n+l)[^/(n+fc 1 +2)(n+fc 1 +l)+^/(n+2)(n+l)] . 



n 2 m 2 



na + ^) 

3=0 



k 2 

na + 

i'=0 



k'2—j' ■ 



(20) 



( 2+ 3)fc 1+ a 



In the framework of the harmonic approximation (i.e. e/n — > where e is an arbitrary finite 
number) the second part of (f20j) reduces to 



ki _ fci-3 _ fc 2 -i 
ui ~ — n 2 m 2 
P 2 



(21) 



Expression (|21|) shows that there are three cases can provide RCP in the evolution of the Qi(T), 
which are: (fci,&2,Mi) = (3,1,3/2) and (/ci,&2,Mi) = (1, 3, 1/2), (2, 2, 1). For the latter cases the 
values of n and fh have to be comparable. Information about these cases has been shown in Figs. 
2(a)-(c) for given values of the parameters. It is obvious that we have three different shapes of 
the RCP. Form Fig. 2(a) revivals and secondary revivals are remarkable having shapes similar to 
those of the (d Z {T)) k 1= k 2 =i (compare to Fig. 1). Fig. 2(b) includes revivals and collapses but 
their forms are different form those of the (<5" z (T))fc 1= fc 2= i. Nevertheless, in Fig. 2(c) the RCP is 
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systematic, i.e. revivals are compact and occur periodically with period tt. Actually, this behaviour 
is quite typical with the evolution of the atomic inversion of the single-mode two-photon JCM, i.e. 
(<r 2 (r))fc 1 =2,fc 2 =o- This can be analytically realized by applying the harmonic approximation to the 
argument of the cosine in ([18]) for the case (fci, k%) = (2, 2) as 

A n+2 , m - K, m = vV + l)(™ + 2) [V(" + 3)(n + 4) - y/{n + l)(n + 2) 

~ ^(m + l)(m + 2)[n + | - n - |] (22) 



= 2^(771 + 1)^ + 2). 

The last line in ([22]) is typical the argument of the cosine of the (<5' 2 (T))fc 1= 2 i fc 2 =o, e.g. see equation 
(12) in [3l| for x = 0. 

Now we deduce the rescaled squeezing factors for the cases (k\, k^) = (3, 1) and (2, 2), which can 
give (<7 2 (T))fc 1= fc 2=1 and (a z (T)) ^=2^2=0^ respectively. From Fig. 2(a) and the discussion given 
above we can write the rescaled squeezing factor for (k\, k^) = (3, 1) as 

Vl(T) ~ — (mW) — • (23) 

Similarly the rescaled squeezing factor for the case (fci, /C2) = (2,2) is 

V , (T) _ <ni(0))-Qi(T) 

Vl{T) - (Mo)) ' ( } 

Expression (|23|) and (|24p have been depicted in Figs. 3(a) and (b) for the values of the interaction 
parameters as those given for Figs. 2(a) and (b), respectively. Comparison between Fig. 1 and 
Fig. 3(a) demonstrates our conclusion. Also Fig. 3(b) provides completely typical shape as that 
of (a z {T)) kl=2M=0 (see Fig. 2 in [32|). 

IV. TWO-MODE SQUEEZING 



In this section, we use procedures similar to those given in section 3 to investigate the RCP in 
the evolution of the two-mode squeezing for TJCM. Starting with the two-mode squeezing factors, 
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which can be expressed as 

S 2 (T) = Re(4(T)ai(T) +ai(T)a 2 (T)) -2Re(ai(T))Re(a 2 (r)) 



2 

+ 

J 



2i 



E [Ma](T)aj(T)) +Re(af(T)) - 2(Re(a J (T))^ 

(25) 

Q 2 (T) = Re(a 2 (T)a!(T) - ai(T)a 2 (T)) - 2Im(a 1 (T))Im(a 2 (T)> 
+ E " Re(af (T)) - 2(lm(d J (T))) 2 " . 

J=1 L 

We discuss the natural and numerical approaches for two- mode squeezing in the following parts. 

A. Natural phenomenon 

Explanations similar to those given in subsection 3.1 the two- mode squeezing factors can give 
direct information on the corresponding (a z (T)) for initial states, which satisfy simultaneously the 
following conditions: 

(a i (T)) = > (a 2 j (T)) = 0, J = 1,2 (26) 

We should stress that (|26p requires the two modes to be initially prepared in such type of states. 
This is different from (|15p of the single-mode squeezing, which requires the mode under consid- 
eration only to be initially in such states. In this case the two squeezing factors (|25p are typical 
and can give (a z (T)). Similar to the single-mode squeezing case this can be verified when the two 
modes are initially prepared in the /-photon coherent states with I = 3,4,.., etc. For such states 
one can easily prove that 

{a z (T)) = (at(O)oi(O)} + (a|(0)a 2 (0)) + 1 - S 2 (T). (27) 



B. Numerical simulation 



Discussion similar to that given in subsection 3.2 one can easily prove that RCP can occur in 
Q 2 (T) when (fei, k 2 ) = (3, 1), (1, 3) and (2, 2). Also one can easily realize for the case (ki, fc 2 ) = (1) 3) 
that Q 2 {T) exhibits RCP, which is the combination from those shown in Figs. 2(a) and (b), i.e. 
RCP is different from that of the (a 2 (r))/ Cl= i j fc 2= i. Nevertheless, for the case (ki,k 2 ) = (2,2), 
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Q2(T) can give information on the (<7 z (T)}fc 1=2> fc 2= o. The rescaled squeezing factor, which is typical 
(<5" 2 (r))fc 1= 2,fc 2 =0) can be evaluated as 

T/ , m (ni(0)) + (n 2 (0))-Q 2 (r) 

V2{T) - (M0)) + (M0)) ■ (28) 

V. SUM AND DIFFERENCE SQUEEZING 

In this section we investigate the occurrence of the RCP in the evolution of the sum and differ- 



ence squeezing factors 33j . In these factors the intermode correlation is involved in the quadrature 
squeezing. We have noted for numerical-simulation approach that the technique given in section 
3 is partially working for sum and difference squeezing. More illustratively, it can give the exact 
values for the transition parameters kj whose squeezing factors exhibit RCP but it fails to pro- 
vide the correct rescaled squeezing factor. Nevertheless, this difficulty can be numerically treated. 
We have noted that the sum and difference squeezing give only information on the occurrence of 
the revivals (not secondary revivals) in (o z (T))k 1= i^ 2 =i- Moreover, sum squeezing can provide 
(<r 2 (r))fc 1= 2 i fc 2 =o, however, difference squeezing fails. We discuss all these results in the following. 



33] 



For the sum squeezing we have 

X = -[aia 2 + a\al], Y = ^[a\a\ — aid 2 ], C = a\ai + a|a 2 + 1. (29) 

Therefore, sum-squeezing factors can be expressed as 

S 3 (T) = Re{a\{T)a 2 2 {T)) + {^(T^T^T)^)) - 2(Re(a 1 (T)a 2 (T))) \ 



(30) 



Q 3 (T) = <aJ(r)oi(r)4(T)oa(r)) - Re(ai(T)a 2 2 (T)) - 2(lm(a 1 (T)a 2 (T))Y 



For difference squeezing the quadratures X and Y can be obtain from those in (|29p by using the 
transformation a\ «-> a\ and consequently C = a 2 a2 — a\a\. For the sake of simplicity we assume 
that the two modes are initially prepared in states having the same photon-number distribution 
with ax = q 2 . This leads to that {a\{T)a 2 {T)) = <a}(T)oi(T)), i.e. \{C{T))\ = 0. Under these 
conditions the difference squeezing factors take the forms 

S A {T) = Re(af (T)a|(T)) + ^(^(T^en^T)) + {a\{T)a x {T)) - 2(Re(d 1 (T)a 2 (T)>) 2 , 

(31) 

Q 4 (T) = {aliT^T)^)^)) + (aJ(r)oi(T)) - Re(a{ 2 (T)a 2 2 (T)) - 2(lm(a 1 (T)a 2 (T))) 2 . 



13 



1.0 n 



0.5 - 




-0.5 - 



— 1.0 i i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i i 

0.0 1.0 2.0 3.0 4.0 5.0 

T 

FIG. 4: The rescaled squeezing factor Va(T) against the scaled time T when the optical cavity modes are 
initially prepared in the three-photon coherent states and (fci, fe, ori, 0-2) = (lj 1> 5, 5). 

A. Natural phenomenon 

Now we seek states that evolve with the TJCM causing the evolution of the (of (T)o|(T)) and 
{a\(T)cL2{T)) close to zero. This can occur if one of the two modes at least is initially in, e.g., the 
three-photon or four-photon states, (cf. ©). For such states expressions (|3U|) of sum squeezing 
reduce to 

Q 3 (T) = S 3 (T) = (aliT^i^aliT^iT)). (32) 
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FIG. 5: The sum-squeezing factor against the scaled time T when the optical cavity modes are initially 
prepared in the coherent states with ct\ = a 2 = 5 for (fa, fa) = (3, 1) (a) and (2, 2) (b). 

In the framework of harmonic approximation the rescaled squeezing factor associated with ()32|) . 
which can provide the corresponding atomic inversion, is 

2<ni(0))(n 2 (0)> + (ni(0)) + (n 2 (0)) + 1 - 2S 3 (T) 



Vs(T) 



(33) 



(ni(0)> + (ni(0)> + 1 

In Fig. 4 we have plotted (j33H when the modes are initially prepared in the three-photon coherent 
states. Actually, we have found that Vs(T) = (a z (T)). From Fig. 4 the revivals and secondary 
revivals are remarkable. Additionally, the revival times of this case are three times smaller than 
those of the initial coherent light since we are dealing with three-photon states (compare Fig. 1 
and Fig. 4). 

Similarly for the difference squeezing (|3ip we can obtain 



((n 1 (0)) + l) 2 -Q 4 (T) 
(ni(0)) + l 



(34) 



B. Numerical simulation 

Similar arguments as those given in subsection 3.2 lead to that if sum-squeezing factor exhibits 
RCP the quantity Re(a|(r)d|(r)) is responsible for this. From (jlUj) and the harmonic approxima- 
tion technique we arrive at 

oo 

{a\{T)a 2 2 {T)) ~ a\a\ ^ C\ m cos[T(A n+2)m+2 - A n , m )]. (35) 

n,m=0 
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Now we seek the proportionality factor, which can be obtained from the following 



2v^r+l)(m + l)' 



After lengthy calculation but straightforward (|36|) reduces to 

* k~i — 1 &2 — 3 fc^ — 1 ^2~5 — 3 &2 — 1 fc^ — 3 &2 — 3 

^2 — 5(2^2^1 2 «2 2 + ^2^1 2 ^2 2 + 2A;in 1 2 n 2 2 + 4kik2 ; n 1 2 n 2 2 



fel — 3 fco — 5 fc^ — 5 &2 — 1 k-^—5 fc2 — 3 fc^ — 5 &2 — 5 

+1k\k\n x 2 n 2 2 + A^r^ 2 n 2 2 + 2k\k2fi x 2 n 2 2 + k\k\n x 2 n 2 2 



(37) 



Similar to the single-mode squeezing case when fi\ ~ Y12 there are three cases, which can provide 
RCP in the evolution of the Qs(T), namely, (k%, fo) = (3, 1), (1, 3) and (2, 2). For these cases the 
proportionality factor is [12 = 2. Actually, this factor cannot give the correct rescaled squeezing 
factor. This fact can be realized from Figs. 5(a) and (b), in which we have plotted the sum 
squeezing factors for (^1,^2) = (3,1) and (^1,^2) = (2,2), respectively. From these figures one 
can see that for the case (k\,k2) = (3,1) sum squeezing can give in principle information on 
fe 1= fc 2 =l! however, for (fci,/^) — (2,2) it can provide the evolution of the {d' z (T)) kl= 2 tk2= o. 
From Figs. 5 and expression (I30p the rescaled squeezing factor is 

= (M0))(n,m-Q 3 (T)_ 
41 ' (n!(0))(n 2 (0)) V ; 

Expression ([38]) gives (a z (T)) kl=k2=1 and (a z (T)) kl=2 ,k 2 =o for (ki,k 2 ) = (3,1) and (2,2), respec- 
tively. 

On the other hand, for the difference squeezing we found that the RCP can be remarked in the 
evolution of the Qa{T) when (ki, fe) = (3, 1), (1, 3), however, for (k\, ^2) = (2, 2) the technique fails. 
To demonstrate these cases we have plotted Figs. 6(a) and (b) for Qi{T) when (fci, fo) = (3, 1) and 
(2,2), respectively. From Fig. 6(a) RCP is established but it is completely different from that of 
the (a z (T)) kl=k2= i, however, Fig. 6(b) exhibits periodically inverted peaks, i.e. it does not exhibit 
RCP. Now from Fig. 6(a) and expression (f3Tj) . the rescaled squeezing factor is 



(ni(0))« ni(0)) + l)-Q 4 (%) 
(«l(0)) 2 



w = (39) 



We have numerically checked (|38|) and ()39[) for the case (ki, k%) = (3, 1) and found that they give 
typical behaviour as that of the {a z (T)) kl=k2= i except that the secondary revivals are absent. 
Additionally, the widths of the revival patterns of (|39p are a little bit greater than those of the 
{cr z {T)) kl=k2=1 . 
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FIG. 6: The difference-squeezing factor against the scaled time T when the optical cavity modes arc initially 
prepared in the coherent states with a\ — ol^ — 5 for (ki, k,2) = (3, 1) (a) and (2, 2) (b). 

VI. CONCLUSIONS 



In this paper we have discussed the possibility of including the squeezing factors of the two- 
mode multiphoton JCM information on the atomic inversion of the standard TJCM. In contrast 
to the single-mode JCM [131 ] we have various types of quadrature squeezing, namely, single-mode, 
two-mode, sum and difference squeezing. Two approaches have been applied for all these types, 
which are natural phenomenon and numerical simulation. Natural approach has been devoted to 
the standard TJCM and found that there is a class of states their squeezing factors provide the 
corresponding atomic inversion. For the numerical-simulation approach we have shown that for 
specific value of the transition parameters, in particular, k\ + ki = 4 the ^-quadrature squeezing 
factor can provide RCP similar to that associated with the (a z (T))i i . 1= k 2= i or (o" z (T))fc 1=2i fc 2 =o 
based on the values of kj. Specifically, for (fci,^) = (3,1) single-mode squeezing factor can 
give information on {a z (T))k 1 =k 2 =ii however, sum and difference squeezing factors give partial 
information. On the other hand, for (fci,^) = (2,2) single-mode, two-mode and sum squeezing 
factors give information on the (a z (T)}k 1= 2,k 2 =o- Also we have deduced the rescaled-squeezing 
factors for all these types giving information on the atomic inversion. 

We conclude this paper by mentioning that the influence of the values of the atomic phases 



on the phenomenon under consideration is the same as that for the single-mode JCM 



In 



other words, for natural (numerical) approach the squeezing factor is sensitive (insensitive) to the 
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initial atomic state, i.e. natural approach can provide "coherent trapping" [34]. Also for numerical- 
simulation approach we have numerically checked that the RCP occurs in the quadrature squeezing 
only for k\ + k<z = 4. 
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